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Modifications of turbulence and turbulent transport associated with a bias-induced 

confinement transition in LAPD 

T.A. Carteri'2'[*]and J.E. Maggs^ 

^Department of Physics and Astronomy, University of California, Los Angeles, CA 90095-1547 
^Center for Multiscale Plasma Dynamics, University of California, Los Angeles, CA 90095-1547 

Azimuthal flow is driven in the edge of the Large Plasma Device (LAPD) [W. Gekelman, et. 
al, Rev. Sci. Instr. 62, 2875 (1991)] through biasing a section of the vacuum vessel relative to 
the plasma source cathode. As the applied bias exceeds a threshold, a transition in radial particle 
confinement is observed, evidenced by a dramatic steepening in the density profile, similar to the L- 
to H-mode transition in toroidal confinement devices. The threshold behavior and dynamic behavior 
of radial transport is related to flow penetration and the degree of spatial overlap between the flow 
shear and density gradient profiles. An investigation of the changes in turbulence and turbulent 
particle transport associated with the confinement transition is presented. Two-dimensional cross- 
correlation measurements show that the spatial coherence of edge turbulence in LAPD changes 
significantly with biasing. The azimuthal correlation in the turbulence increases dramatically, while 
the radial correlation length is little altered. Turbulent amplitude is reduced at the transition, 
particularly in electric field fluctuations, but the dominant change observed is in the cross-phase 
between density and electric field fluctuations. The changes in cross-phase lead to a suppression 
and then apparent reversal of turbulent particle flux as the threshold is exceeded. 
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PACS numbers: 52.35.Ra, 52.25.Fi, 52.35.Kt 



INTRODUCTION 



Turbulence has long been recognized as the dominant 
cause of cross-field transport in magnetically-confined 
plasmas [T Significant progress in the effort to confine 
plasmas for fusion energy came with the discovery of the 
high- confinement mode or "H-mode," where the sponta- 
neous formation of an edge transport barrier occurs after 
the heating power exceeds a threshold [2 . A key char- 
acteristic of these transport barriers is the presence of a 
localized, non-uniform radial electric field with associated 
cross-field E x B flow and flow shear [3l |4] . The impor- 
tant role of radial electric field in the confinement transi- 
tion was first demonstrated clearly in experiments in the 
Continuous Current Tokamak (CCT) O |6] where bias- 
ing was used to establish a radial electric field and trigger 
an H-mode confinement transition. Since these pioneer- 
ing experiments, biasing has been used on a number of 
confinement devices to study transport barrier formation 
and turbulence modification [7, 8 , l9l fTO) ITT| 1121 [T3] . 

Flow and flow shear can modify turbulence and trans- 
port through a number of processes. First, the shear 
can lead to radial decorrelation or "shearing apart" of 
turbulent eddies, reducing their radial extent and there- 
fore transport effectiveness [14 . Also, the amplitude of 
turbulent fluctuations can be reduced through the mod- 
ification of the nonlinearly saturated level of the turbu- 
lence |3 . Finally, electrostatically-driven particle trans- 
port is dependent on the correlation between density 
and electric field fluctuations in low frequency turbulence 



(where E x B is the dominant particle response): 
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where n{u;) and Eoiuo) are the Fourier transforms of 
density and azimuthal electric field, ^n^Ee is the cross- 
coherency between n and Eq^ and On^Ee is the cross- 
phase between density and electric field fluctuations [15] . 
Flow and flow shear can modify both the cross-phase 
and cross-coherency between density and potential fluc- 
tuations, leading to a reduction of transport even in the 
absence of amplitude reduction [161 llZ] • 

Changes in turbulence and transport associated with 
azimuthal flow have been investigated in the Large 
Plasma Device (LAPD) at UCLA. Flows are driven 
through biasing the vacuum chamber wall relative to the 
plasma source [TH^, |19[ |20]. Excitation of the Kelvin- 
Helmholtz instability has been observed in LAPD when 
a strong flow gradient is imposed (through modification 
of the magnetic topology) [191 IIO]- For flow driven in 
a straight magnetic field line configuration in LAPD (re- 
sulting in smaller but non- negligible flow shear) , a transi- 
tion in particle confinement is observed. Models of diffu- 
sive particle transport and changes in transport are pre- 
sented in an earlier publication [18 where it is reported 
that transport rates change from Bohm to classical due 
to edge rotation. The observation of a confinement tran- 
sition in a simple magnetic geometry provides the op- 
portunity to perform a detailed study of modifications of 
turbulence and turbulent transport in a system free from 
complications associated with toroidal systems (e.g. field 



line curvature, poloidal asymmetry, trapped particles). 
Good diagnostic access to LAPD provides for detailed 
measurements of the spatial and temporal characteristics 
of the turbulence. 

Here we summarize the primary results reported in this 
paper. Measured turbulent transport flux is reduced and 
then suppressed, leading to a confinement transition, as 
the applied bias is increased. The threshold in the ap- 
plied bias is linked to radial penetration of the driven 
azimuthal fiow. Two-dimensional measurements of the 
turbulent correlation function show that the azimuthal 
correlation increases dramatically during biasing, with 
the high m-number modes involved in turbulent trans- 
port becoming spatially coherent. However no significant 
change in the radial correlation length is observed asso- 
ciated with the confinement transition. As the bias is in- 
creased above threshold, there is an apparent reversal in 
the particle fiux (indicating inward transport). The peak 
amplitude of density and electric field fiuctuations do de- 
crease, but the reduction is only slight and does not fully 
explain the transport fiux reduction. The cross-phase 
between density and electric field fluctuations changes 
significantly as the threshold for confinement transition 
is reached, and explains the reduction and reversal of the 
measured transport fiux. The dynamics of the transition 
have also been studied and show a correlation between 
transport suppression and the overlap of the fiow (shear) 
profile and the density gradient profile. 



II. EXPERIMENTAL SETUP 

The experiments were performed in the upgraded 
Large Plasma Device (LAPD) [21], which is part of the 
Basic Plasma Science Facility (BaPSF) at UCLA. The 
vacuum chamber of LAPD is 18 m long and 1 m in di- 
ameter and surrounded by solenoidal magnetic field coils. 
The plasma is generated by a cathode discharge [22 . The 
cathode is 73 cm in diameter and is located at one end of 
the vacuum chamber, as shown schematically in Fig. [T] 
A molybdenum mesh anode is situated 50 cm away from 
the cathode. A bias of 40-60V is applied between the 
cathode and anode using a solid-state switch [23], result- 
ing in 3-6kA of discharge current. An important aspect 
of the operation of the plasma source is the generation of 
primary electrons with energy comparable to the anode- 
cathode bias voltage. The mesh anode is 50% transpar- 
ent, allowing half of these primaries to travel down the 
magnetic field into the main chamber leading to ioniza- 
tion and heating of the bulk plasma. The cathode is 
situated at the mouth of the solenoid and therefore sits 
in a region of flaring magnetic field. The fiared magnetic 
field maps the 73 cm diameter cathode to a ~ 56 cm di- 
ameter region in the main chamber. Primary electrons 
from the source are isolated to this region in the cham- 
ber and lead to a fast electron tail for r < 28cm. Typical 
plasma parameters in LAPD discharges are ng ^ 5 x 10^^ 
cm-^, Te - 7 eV, T^ - 1 eV, and B < 2kG. In the 



experiments reported here, the primary plasma species 
was singly-ionized helium and the magnetic field strength 
used was 400G. Modeling [18 and spectroscopic measure- 
ments show that the ionization fraction of the LAPD 
plasma is > 50% and therefore Coulomb collisions are 
the most important collisional process. However, neutral 
collisions (charge exchange in particular) are important 
for establishing the radial current during biasing. 



Discharge Circuit 




FIG. 1: Schematic of LAPD including wall biasing circuit. 

Measurements of density, temperature, fioating po- 
tential and their fiuctuations are made using Langmuir 
probes. A four-tip probe with 0.76mm diameter tanta- 
lum tips arranged in a diamond pattern is used as a triple 
Langmuir probe and particle fiux probe. Two tips are 
separated 3mm along the field and are used as a double 
probe to measure ion saturation current (/gat oc ne^/T~e). 
The remaining two tips are separated 3mm perpendicular 
to the field and measure fioating potential for deriving az- 
imuthal electric field fiuctuations and temperature using 
the triple Langmuir probe method [24 . Radial parti- 
cle transport can be evaluated directly using measured 
density and electric field fiuctuations through Eqn. [l] 
Flows are measured using a Gundestrup (Mach) probe 
with six faces [25] . The fiow measurements are corrected 
for the finite acceptance angle of the probe faces, which 
are smaller than the ion gyroradius [26 . 

Using the difference in fioating potential to determine 
the azimuthal electric field is problematic in the pres- 
ence of fast electron tails, such as exist on field lines con- 
nected to the cathode source. In the presence of fast 
electrons, differences in fioating potential may no longer 
be proportional to differences in plasma potential and 
thus may not accurately measure the electric field. In 
addition, the fast electron tail, which sets fioating po- 
tential, is much less collisional than the bulk plasma so 
that the fioating potential is no longer a localized quan- 
tity. This de-localization skews the correlation with lo- 
cally measured density fluctuations and thereby affects 
the flux measurement. Measurements of /gat are made 
with a double-Langmuir probe biased to 70V to reject 
primary electrons, and therefore should not be affected 
by the fast electron tail. Thus, in this study, we report 
on the properties of ion saturation fluctuations every- 
where in the plasma column. Properties of electric fleld 



fluctuations and cross-correlation flux measurements are 
presented everywhere in the plasma column, but a cau- 
tion is issued in regard to measurements made on fleld 
lines where primary electrons are present (r < 28cm). 
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FIG. 2: (color online) Radial profiles of (a) floating potential 
and (b) perpendicular Mach number for no bias and 175V 
bias. Note the radial location of field lines that connect to 
the cathode edge are marked with "CE" in this figure and 
subsequent figures. 

The edge plasma in LAPD is rotated through bias- 
ing the vacuum vessel wall positively with respect to the 
source cathode. Figure [l] shows a schematic of the LAPD 
along with the rotation bias circuit, which includes a 
capacitor bank and an IGBT (Insulated Gate Bipolar 
Transistor) switch. The rotation bias is pulsed, with the 
switch closed for 4 ms during the steady portion (current 
flattop) of the LAPD discharge. Typical bias voltages are 
on the order of lOOV, and result in around lOOA of peak 
radial current. The radial current provides the torque 
(through J X B forces) needed to spin up the plasma. 
The measured current is consistent with ions carrying 
the radial current through ion-neutral collisions (Peder- 
son conductivity) [iSj . The ions carry the current radially 
until they reach fleld lines which connect to the cathode 
in the plasma source, where the current closes via par- 
allel electron currents. The influence of the applied wall 
bias is therefore limited to regions of the plasma that 
are not connected along the magnetic fleld to the cath- 
ode source. This is demonstrated in Figure [2] (a), which 
shows the measured floating potential as a function of ra- 
dius for no bias and a bias of 175V. Measurements shown 
in this flgure (and subsequent flgures) are made with a 
1 cm radial resolution through moving a single probe on 
a shot-to-shot basis. The trace shown is the result of 
a time average over 1.3 ms of data during the flat-top 
of the LAPD discharge current pulse. Field lines in the 
region r < 28 cm connect to the cathode, and the float- 
ing potential is unchanged by the applied bias in this 



region. Note that in this flgure (and many following flg- 
ures) the radial location of the fleld lines which connect 
to the cathode edge is marked with "CE." Measurements 
of azimuthal flow (Mach number) for the same case are 
shown in Fig[2](b), indicating that driven flows are also 
localized to fleld lines not connected to the plasma source. 



III. EXPERIMENTAL RESULTS 
A. Threshold for confinement transition 

A radial particle conflnement transition is observed 
with biasing [18 . Figure [s] shows proflles of plasma 
density, ion saturation current (/gat) fluctuations, and 
perpendicular Mach number for two cases: unbiased 
and 175V bias. Without bias, radial turbulent parti- 
cle transport causes the density proflle to extend well 
past the cathode edge, with a fairly gentle gradient 
{Ln = |Vlnn|~^ ~ 10 cm), as shown in Fig. p](a). For 
comparison, the ion sound gyroradius is ps ~ 1.2 cm (he- 
lium, Tg ~ 6eV, B = 400G) and the ion gyroradius is 
pi ~ 0.5 cm {Ti ~ leV). In the biased case, the mea- 
sured density proflle steepens dramatically (L^ ~ 4 cm) . 
Detailed transport modeling of the LAPD plasma shows 
that the proflle before biasing is consistent with radial 
diffusion at the Bohm rate. The steepened proflle after 
biasing is consistent with classical collisional diffusion, a 
factor of ~ 100 drop in the diffusion coefflcient [18 . In 
the rotating plasma, the density proflle is essentially set 
by the plasma production source proflle since the radial 
transport rates are so slow. 

The electron temperature proflles measured using the 
triple probe technique are shown in Fig. 3 (b). The 
electron temperature drops rapidly past the edge of the 
cathode. The electron beam from the cathode acts pri- 
marily as a heat source once the plasma is formed. The 
fully formed plasma column is hot and dense enough 
that the tail of the Maxwellian distribution dominates 
the beam electrons in the ionization process. In the ro- 
tating plasma the electron temperature is elevated in the 
region outside the plasma source, perhaps due to fric- 
tional heating from neutral drag. The Mach-probe mea- 
sured azimuthal flow proflle, with and without bias, is 
repeated in Fig. [s] (c) for reference. 

The observed proflle steepening occurs for bias volt- 
ages above a threshold value. Figure [4] shows the mini- 
mum density gradient scale length (L^) versus bias volt- 
age. For biases above ~ lOOV, the density gradient 
scale length decreases dramatically, indicating a transi- 
tion in the radial particle conflnement. As the bias is in- 
creased, the density gradient scale length saturates (near 
Ln ^ ^ cm), consistent with modeling which shows that 
transport beyond the conflnement transition threshold is 
suppressed to the irreducible lowest level set by classical 
transport [18]. 

The threshold for conflnement transition in LAPD ap- 
pears to be set by radial penetration of the bias-driven 



3.0 
2.5 
2.0 
1.5 
1.0 

0.5 
0.0 
10 

8 



^. 4 




1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 



f (a) 


— No Bios 

- - ■ 1 75V Bios 


-i 


"-^l- - 






ce;^"^-^ 


-i 


^ (b) 




r (^^ 




J 


/^ 


-i 


-;r<^ — 


CEI " ■ — 


■ 


. 1 . : 


-: 



15 20 25 30 35 40 45 

r (cm) 

FIG. 3: (color online) Profile of (a) density, (b) Electron 
temperature, and (c) Mach number before and after a con- 
finement transition. 
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FIG. 4: Minimum radial density gradient scale length as a 
function of applied bias. (This figure reprinted from refer- 
ence ^8j) 



flow. Figure [sja) shows the radial profile of azimuthal 
flow for several values of the bias, including a fit to 
the data of the form a\ tanh((r — a^)!^^) + <^4^^ + <^5- 
Figure [5|b) shows the profile of the shearing rate, 7 = 
dvo/dr^ derived from the fit. This technique is used to 
reduce the noise introduced by taking radial derivatives 
of the raw data. Below the transition threshold (75V) 
the flow and flow shear is concentrated in the far edge 
plasma, and has not penetrated in to the cathode edge 
at r ~ 28cm. As the threshold is reached (lOOV) the 
peak value of the flow shear has not increased beyond 
the below-threshold value, but the flow has moved ra- 



dially inward so that some finite flow shear exists near 
the cathode (plasma source) edge. As the bias is fur- 
ther increased, the flow amplitude and shear increase, 
but the shape of the profile is relatively unchanged. The 
observed transition is therefore not due to systematic in- 
crease of the shearing rate through a threshold value, but 
is instead associated with a sudden appearance of a large 
shearing rate at the plasma source edge. 
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FIG. 5: (color online) For several bias values: (a) Radial 
profile of azimuthal velocity, with fit. (b) Radial profile of 
shearing rate, derived from fit to azimuthal velocity profile. 

If the change in radial transport rates is associated 
with shear suppression of turbulent transport, the shear- 
ing rate should be large enough to compete with char- 
acteristic timescales of the turbulence, often taken to be 
the linear growth rate or the nonlinear eddy turnover 
time. The turbulence in the edge of LAPD is known 
to be drift- Alfven in nature [27 , and the characteris- 
tic frequency of drift waves is the diamagnetic drift fre- 
quency, uod = kyVth,epe/Ln (whcrc ky is the azimuthal 
wavenumber, Vth,e is the electron thermal speed and p^ 
is the electron gyroradius). The drift frequency will be 
taken as a proxy for the linear growth rate, although 
the growth rate of resistively-driven modes is often only 
a small fraction of this frequency [28 . The nonlinear 
timescale in the turbulence, the eddy turnover time, can 
be estimated from measurements of the autocorrelation 
time in the turbulence, Tac- Figure [6] shows the radial 
profile of shearing rate, r~^^ and ujd for three bias val- 
ues: 75V (below threshold), lOOV (near threshold) and 
175V (above threshold). In the calculation of cojj^ the 
effective azimuthal wave number of the fluctuations, ky, 



is taken from spatial correlation measurements (detailed 
later in this paper) which show a dominant azimuthal 
mode number of m ~ 10. The autocorrelation time is 
derived from the lab-frame autocorrelation function com- 
puted from ion saturation current data, which is fit to 
the function ao cos (air) exp (— a2 |r|), where Tac = 1/^2 • 
At the threshold for transition, the shearing rate in the 
edge region is comparable to the local diamagnetic drift 
frequency and the autocorrelation frequency (1/rac). In 
a flowing plasma, the autocorrelation time in the lab 
frame is affected both by flow of turbulent structures past 
the probe as well as temporal decorrelation of turbulent 
structures in the plasma frame (the eddy lifetime) |29 . 
The autocorrelation time due to flow of structures past 
the probe can be estimated as Ty ~ '^^x/'^Aow, where wj_ 
is the azimuthal scale size of the turbulence (azimuthal 
correlation length) and '^aqw is the azimuthal flow veloc- 
ity. An estimate of t~^ is shown for a bias of 175V in 
Fig. IgFc), using an azimuthal scale size of w± ~ 10cm 
(motivated by spatial correlation measurements). The 
r~^ estimate is comparable to the calculated r^^^, indi- 
cating that flow of structures past the probe is a major 
contributer to establishing the lab-frame autocorrelation 
time, and that it can not be interpreted as a measure of 
the eddy turnover time. However, Tac does set a lower 
limit for the eddy turnover time and the actual plasma- 
frame eddy lifetime must be longer than Tac (and there- 
fore r~^. must be smaller than r~^^). It is therefore rea- 
sonable to argue that the flow shear is large enough to 
be dynamically important and thus could suppress turbu- 
lent transport, leading to the observed confinement tran- 
sition. 
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B. Changes to spectrum and transport flux under 
biasing 

Detailed measurements of turbulent spectra and tur- 
bulent particle flux have been performed to evaluate 
changes associated with the confinement transition. Fig- 
ure [7| shows the radial profile of rms /gat fluctuations, 
contours of the /gat power spectrum, and the profile of 
perpendicular Mach number for several values of the ro- 
tation bias. All are time-averaged over the length of 
the bias pulse. The normalized fluctuation amplitude 
is SIsat/Isat ^ 33% at the peak of the fluctuation am- 
plitude in the unbiased case. The fluctuation spectrum 
is broadband prior to biasing, with no evidence for co- 
herent modes. Below threshold (75V), some reduction 
of fluctuation amplitude is seen in the far edge, as well 
as a Doppler upshift in the fluctuation spectrum (espe- 
cially associated with the strong flow feature near the 
wall at r = 50cm). As the transition threshold is ap- 
proached (lOOV) and exceeded (150V), the fluctuation 
amplitude is reduced and concentrates on the steepened 
density gradient, but the spectrum remains broadband. 
The rearrangement of the fluctuation profile accompa- 
nies changes in the plasma density profile (not shown). 



FIG. 6: (color online) Radial profiles of the shearing rate 
(Texb), the turbulent autocorrelation time (rac), and the drift 
frequency (cod) for three bias values: 75V (below threshold), 
lOOV (threshold) and 175V (above threshold). 



A strong Doppler shift is also observed in the spectrum 
localized to the region of strong azimuthal flow. 

Radial decorrelation, or reduction of the radial size 
of turbulent eddies, is one possible mechanism by which 
turbulent transport might be suppressed by sheared flow. 
Measurements of the two-dimensional (cross-field) turbu- 
lent correlation function have been performed to evaluate 
the role of this mechanism in the confinement transition. 
Figure |8] shows measurements of the two-dimensional cor- 
relation function of fluctuations in ion saturation current 
for unbiased and strongly biased (220V) cases. These 
measurements are made using two probes separated along 
the magnetic field by 0.69m. One probe is fixed in space 
near the cathode edge (at r = 26cm) and used as a refer- 
ence to correlate against, while the second probe is moved 
shot-to-shot to make measurements at 775 locations in 
the plane perpendicular to the background magnetic field 
(the density gradient is in the —x direction). Prior to bi- 
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FIG. 7: (color online) Radial profile of RMS /sat fluctuations, 
Contours of /sat power spectrum versus radius and frequency, 
and Mach number profile for several bias values. The radial 
profiles of RMS fluctuation level and Mach number in the 
unbiased case are shown as dashed lines for comparison with 
the biased cases. 



asing, the correlation function is nearly isotropic, with a 
slightly longer correlation length in the radial direction 
{x in the figure). The correlation lengths in the unbi- 
ased case (radial and azimuthal) are < 5 cm, around 
5 times the ion sound gyroradius and 10 times the ion 
gyroradius. Under biasing, the azimuthal {y) correlation 
length increases dramatically. However, the radial extent 
of the correlation function is not reduced significantly 
and therefore no strong evidence for radial decorrelation 
is found. In addition, a clear pattern of maxima and 
minima are visible in the correlation pattern measured 
in the rotating plasma. This pattern resembles the pat- 
tern expected of a coherent, high m- number cylindrical 
eigenmode. 

There are two potential explanations for transport 
reduction based on the measured correlation function. 
First, because the correlation function is strongly elon- 
gated in the presence of the flow, the eddy turnover time 
is lengthened: it takes longer for a blob of plasma to 
travel from one side of the eddy to the other in the ra- 
dial direction. So the effective transport time step is 
lengthened {i.e. the eddy turnover time increases), even 
though the step size (the radial correlation length) is 
largely unchanged. This explanation seems inconsistent 
with Fig. [6J which shows the autocorrelation time de- 
creasing with increasing bias. However, the lab-frame 
autocorrelation time may be dominated by flow effects. 
Second, the spatial coherency of the fluctuations has in- 
creased dramatically as a result of the biasing. Turbu- 
lent particle diffusion is effective only in the presence of 
decorrelation: coherent vortices can flatten the density 
profile over their width through mixing, but do not lead 
to radial diffusion. 

The expression for particle flux given in Eq. [l] can be 
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FIG. 8: (color) Two-dimensional cross- field correlation func- 
tion for /sat fluctuations in unbiased and biased plasmas 



directly evaluated using data obtained by a particle flux 
probe. Ion saturation current is measured to determine 
density fluctuations while floating potential is measured 
at two azimuthally separated probe tips in order to esti- 
mate fluctuations in the azimuthal electric field Eq. Fig- 
ure [9] (a) shows the radial profile of turbulent particle 
flux measured for four different bias voltages. The flux 
is computed using the FFT of ion saturation current and 
electric field fluctuations, taking the final 2.6 ms of the 
4 ms bias pulse as the window for the FFT computa- 
tion. The measured particle flux is reduced and then 
suppressed as the bias is increased to llOV (just past the 
threshold voltage). As the bias is further increased, there 
is a reversal in the measured particle flux indicating in- 
ward transport. The particle flux in Eq. [l]is established 
by the fluctuation amplitudes and the cross-coherency 
and cross-phase between density and electric field fluctu- 
ations. Figure [9] (b) and (c) show the radial profiles of 
density and electric field fluctuations for the same bias 
voltages. The density fluctuations become more localized 
in space (peaked on the steepened density gradient), but 
the peak amplitude does not decrease significantly. The 
electric field fluctuations in the gradient region are also 
slightly reduced after the threshold is exceeded. How- 
ever, electric field fluctuations grow stronger in the far 
edge with increased bias. As noted earlier, measurements 
of electric field fluctuations (and hence computed parti- 
cle flux) are suspect in the presence of primary electrons 
(r < 28cm). The radial region in which primary electrons 
are present is specifically marked in Fig. |9](a) and (c). 

The fast Fourier transform (FFT) power spectrum of 
/sat fluctuations is shown in Fig. [lOja) (averaged over 
spatial locations 28cm < r < SOcmyT The spectrum is 
broadband before and after the transition, with decreased 
power at low frequency and a Doppler upshift in the 
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FIG. 9: (color online) (a) Particle flux profile, (b) root-mean- 
square (RMS) /sat fluctuation profile, and (c) RMS Eg fluc- 
tuation profile for several bias values. 



spectrum apparent with increasing rotation. The cross- 
coherency between density and electric field fluctuations 
actually increases with bias, as shown in Fig. [To](c). Al- 
though the electric field fluctuation amplitude does de- 
crease, the dominant cause of the observed flux behav- 
ior with bias appears to be due to the cross-phase. As 
shown in Fig. 10 (d), the cosine of the cross-phase is very 
favorable for outward transport below the threshold, is 
modified to near zero around the threshold, and finally 
reverses sign at larger bias values. 

Figure [TT] shows a comparison between the radial pro- 
files of the measured flux and the flux determined from 
transport analysis ^. In the unbiased, non-rotating 
case, the measured and predicted fiux agree remarkably 
well on field lines that do not connect to the cathode 
(r > 28cm). However, on cathode connected field lines 
(where primary electrons exist), the measured flux is near 
zero. One possible explanation for these observations is 
the influence of 50eV primary electrons from the source 
on the floating potential measurement. The fastest elec- 
trons establish floating potential, and on the cathode- 
connected field lines these are the primary electrons and 
not the bulk electrons (Tg ~ 5eV). The primaries may not 
participate in the drift- wave dynamics or if they do, they 




f (kHz) 

FIG. 10: (color online) For the same bias values, averaged 
over 28 < r < 30cm and using a 2.6 ms time window at 
the end of the bias pulse for the FFT: (a) /sat FFT power 
spectrum, (b) Eg FFT power spectrum, (c) Cross-coherency 
and (d) cross-phase between /sat and Eq. 



have extremely long mean free path (comparable to the 
machine size) and any perturbation imparted on them 
might be phase-mixed away as they are highly nonlocal 
in the wave. Figure [9] is consistent with this argument, 
showing that the electric field fluctuations peak outside 
of the cathode radius while fluctuations in /gat (which is 
measured with Vbias ^ 70V to reject primary electrons) 
peak further inward. 

Inward flux (up the density gradient) is measured in 
these experiments at larger values of wall bias, while the 
transport model flux of the rotating plasma is nearly zero 
(classical transport). Inward flux has been reported in 
other biasing experiments on toroidal confinement de- 
vices [121 Eq] and also in the presence of turbulently 
driven large-scale shear flows in a cylindrical device J3T] . 
In addition, turbulence-driven particle pinches are pre- 
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FIG. 11: Comparison between model and measured fluxes. 
Primary electrons may influence the ability to measure trans- 
port flux, and the region where they exist is marked in the 
figure. 



dieted theoretieally in the presence of strongly sheared 
flow [32 . However, if the turbulence is driven purely 
by density gradients, only local reversal of the flux is 
allowed, and not reversal across the whole profile, as ob- 
served here (thermodynamically, the density gradient can 
not steepen itself) . The growth of electric field fluctua- 
tions in the far edge (away from density gradients, see 
Fig. [9]) at large bias may be due to a new instability 
driven by other free energy sources such as the flow. It 
is possible that these new fluctuations could drive the 
observed inward flux. However, the measured inward 
flux can not be consistent with transport modeling un- 
less a source of plasma is present in the far edge region 
(r > 28cm) to feed the measured large inward flux. No 
primary electrons are present in this region, but ioniza- 
tion associated with rotation-heated bulk electrons might 
provide an ionization source. The source of the apparent 
fiux reversal in these measurements will be the subject 
of future work. 



C. Dynamics of the confinement transition 

The temporal behavior of the azimuthal flow, turbu- 
lence, and flow profile shape during the confinement tran- 
sition is interesting and provides further support for the 
conclusion that the radial flow profile is central to the 
transport reduction. In this section we focus on the rela- 
tion of the shear rate profile to the density gradient pro- 



rms fluctuation amplitude (/gat fluctuations), and shear- 
ing rate versus radius and time for a case where a con- 
flnement transition occurs (125V bias). The 125V bias is 
applied from t = 7.5ms to t = 11.5ms (both marked with 
dotted lines in the flgure). After the bias is turned on, 
the flow shear penetrates into the cathode edge, where 
it peaks strongly at t ~ 8ms. The flrst effect of the flow 
penetration is an immediate reduction in the fluctuation 
amplitude in the edge region (30 < r < 35). The negative 
correlation between flow shear and edge fluctuation am- 
plitude is especially apparent near the flrst strong peak 
in the shear {t ~ 8ms). The fluctuations further into 
the core plasma are suppressed more slowly, but reach 
their lowest levels 1.5ms after the bias is turned on. On 
a similar timescale, the steepening of the edge density 
proflle occurs and the core proflle flattens. It is therefore 
likely that the reduction of core fluctuations is at least in 
part due to changes in the linear drive of the drift- Alfven 
waves due to core density proflle flattening. The steep- 
ening of the edge density gradient occurs on a timescale 
consistent with the transport timescale set by end losses 
along the straight magnetic fleld [18]. 
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flle. Figure 12 presents contour plots of density gradient. 
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FIG. 12: (color online) Contour plots of density gradient, rms 
/sat fluctuation amplitude and shearing rate versus radius and 
time for 125V bias. 

After the initial transient, the density proflle steep- 
ens and the shearing rate settles into a proflle which 
persists in a quasi-steady state (from t ^ 9ms to t ^ 
10.5ms). The fluctuation amplitude is reduced from the 
pre-transition value, and is localized to the edge gradient 
region. However, late in the bias pulse (near t ~ 10.5ms) 
the fluctuation amplitude grows up again and the density 
gradient reduces. This observation can be explained by 
considering the overlap of the shear and density gradient 
proflles. From Fig. 12 (b) and (c) it is apparent that late 
in the bias pulse the location of peak shear is slowly mov- 
ing radially outward while the location of peak density 
gradient is moving radially inward. This is made clearer 



in Fig. 13 which shows the radial profile of density gra- 
dient and shearing rate, derived from fits, at two points 
in time. To emphasize profile shape and overlap, the 
profiles are normalized to their peak values. Early in the 
quasi-steady phase (t =8. 9ms), significant overlap is seen 
between the density gradient and shearing rate profiles. 
However, later in the bias pulse (t = 10.5ms) the two 
profiles have moved apart and the amount of overlap is 
reduced substantially. Because of the reduced overlap be- 
tween the density gradient (where linear drive for the tur- 
bulence exists) and the shearing rate, shear suppression 
of the turbulence should be less effective. This reduced 
suppression ability can explain the return of fluctuations 
and turbulent transport with the resultant reduction in 
the density gradient. Note that, since similar arguments 
could be made for the overlap of the density and flow pro- 
files (these are the radial integrals of the profiles shown) , 
it is not clear which profiles are determinant. However, 
it certainly appears that the degree of overlap between 
the profiles, whether density gradient and shear rate or 
density and flow, is crucial to the stabilization process. 
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FIG. 13: (color online) Profiles of density gradient and shear- 
ing rate, derived from fits to the density and flow, early and 
late in the bias pulse. 



IV. CONCLUSIONS 

Biasing the central chamber wall of the LAPD, pos- 
itive with respect to the cathode, produces a rotation 
layer on field lines external to the cathode. This rotation 
layer leads to a significant reduction in radial transport 
of plasma. The reduction is accompanied by changes in 
the radial temperature and density profiles of the plasma 
column. In particular the density gradient scale length 
is greatly reduced near the source edge due to the re- 
duction in radial transport rates. For a fixed bias pulse 
length, a threshold for reduction in radial transport ex- 
ists. The threshold voltage is apparently determined by 
radial penetration of the flow into the gradient region of 
the plasma. Once the flow penetrates, a large flow and 
flow shear are immediately present in the gradient region 
and therefore a threshold value for the shearing rate can 
not be determined from this experiment. 

Without rotation, ion saturation current fluctuations 
are broadly distributed from the core to the plasma wall 
with peak amplitude very near, but inside, the cathode 
edge. As plasma rotation proceeds, the fluctuation profile 
adjusts with the changing density profile of the plasma. 
The fluctuations become concentrated in the gradient re- 
gion of the plasma and the peak amplitude is reduced 
by a few tens of percent from the non-rotating value. 
The most dramatic change in the ion saturation fluctu- 
ations is noticed in the cross-fleld correlation function. 
In the non-rotating plasma, the cross-fleld correlation 
function is very nearly isotropic. In the rotating plasma, 
the radial correlation length is about the same as in the 
non-rotating plasma. However, the azimuthal correlation 
length of the rotating plasma is much longer than the ra- 
dial correlation length (by about a factor of flve), and 
the correlation function takes on a noticeable structure, 
as might be expected for a coherent cylindrical eigenfunc- 
tion. 

The most signiflcant change in the cross-correlation be- 
tween azimuthal electric fleld fluctuations (floating po- 
tential) and density fluctuations (ion saturation current) 
occurs in the cross-phase in the frequency range below 
10 kHz (about a tenth of the ion gyrofrequency ) . The 
cross-phase changes systematically with increasing rota- 
tion bias. In the unbiased, non-rotating plasma the cross- 
phase is near unity at low frequencies. As the rotation 
bias increases the cross-phase decreases toward zero and 
then progresses towards minus one for the fully rotating 
plasma. The measured radial fluxes calculated from the 
cross-correlation given in Eq.[l] agree with previous mod- 
eling results for the non-rotating plasma on fleld lines 
outside the cathode, i.e., in the absence of fast tails. In 
the rotating plasma the cross-correlation technique pre- 
dicts negative (inward) fluxes that are inconsistent with 
modeling. 

The change in radial transport depends upon the de- 
gree of spatial overlap between the shear rate and density 
gradient proflles as demonstrated by the behavior near 
threshold. The recovery of radial transport and subse- 
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quent increase in the scale length of the density profile, 
as observed late in the 125V bias case, is accompanied 
by a decrease in the overlap between these two profiles. 
Thus the stabilization process depends upon the spatial 
distribution of the azimuthal flow. 

In these experiments, azimuthal flow was driven exter- 
nally, and the interaction between the driven flow and 
spontaneously driven turbulence was studied. A criti- 
cally important, unresolved issue in confinement tran- 
sitions (e.g. L- to H-mode transition) is the mecha- 
nism behind the spontaneous generation of radial elec- 
tric field and associated sheared flows. Many mech- 
anisms have been proposed, including non-ambipolar 
edge transport, neoclassical effects, and turbulent pro- 
cesses [41 |33l [Ml EH El [37]. In particular, the self- 
generation of flows by turbulence (e.g. zonal flows) has 
recieved a great deal of attention as a likely candidate 
to explain spontaneous confinement transitions [4, 37j. 
Direct evidence for Reynolds-stress-driven flow has been 
reported in a cylindrical device [38] and evidence for en- 
hanced mode-coupling prior to the L- to H-transition 
has been seen in DIII-D [36l [39] (however, similar. 



but not exhaustive, measurements in NSTX did not re- 
veal enhanced mode-coupling [40 ). Spontaneous az- 
imuthal flows are observed in LAPD prior to biasing (see 
Fig. [2|b)) and future work will explore the role of in- 
verse cascade and mode-mode coupling on generation of 
flows and flow-turbulence interaction in LAPD. Future 
work will also consider the role of sheath-driven processes 
present in open-magnetic-field-line configurations in flow 
generation J4T] . 
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